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A Chern-Simons coupling of a new scalar field to electromagnetism may give rise to cosmological 
birefringence, a rotation of the linear polarization of electromagnetic waves as they propagate over 
cosmological distances. Prior work has sought this rotation, assuming the rotation angle to be 
uniform across the sky, by looking for the parity- violating TB and EB correlations a uniform rotation 
produces in the CMB temperature/polarization. However, if the scalar field that gives rise to 
cosmological birefringence has spatial fluctuations, then the rotation angle may vary across the 
sky. Here we search for direction-dependent cosmological birefringence in the WMAP-7 data. We 
report the first CMB constraint on the rotation-angle power spectrum C£ a for multipoles between 
L = and L = 512. We also obtain a 68% confidence- level upper limit of y/Cg**/ (47r) < 1° on the 
quadrupole of a scale- invariant rot at ion- angle power spectrum. 



I. INTRODUCTION 

In this work, we use the cosmic microwave back- 
ground (CMB) temperature and polarization maps of the 
Wilkinson Microwave Anisotropy Probe (WMAP) 7-year 
data release [1 to search for direction-dependent cos- 
mological birefringence (CB). CB is a postulated rota- 
tion of the linear polarization of photons that propagate 
through cosmological distances [2]. It is present, for ex- 
ample, in models where a Nambu-Goldstone boson plays 
the role of quintessence [3], but also in models with new 
scalar degrees of freedom that have nothing to do with 
quintessence [4-7 j. The rotation of the polarization is a 
consequence of the coupling of a scalar field to the elec- 
tromagnetic Chern-Simons term, such that the rotation 
angle a is proportional to the total change A(j) of the field 
(j) along the photon's path. 

Prior to this work, a rotation angle a that is uni- 
form across the sky had been sought in the CMB [8], 
where it would induce parity- violating TB and EB tem- 
perature/polarization correlations [9]. CB has also been 
sought in quasar data [2j [10] . The tightest constraint cur- 
rently comes from a combined analysis of the WMAP, 
Bicep [TT] , and QUAD experiments [12]; it is —1.4° < 
a < 0.9° at the 95% confidence level [T3] . 

There are, however, a number of reasons to expand the 
search and look for a CB angle a(n) that varies as a func- 
tion of position n on the sky. To begin with, a dynamical 
field (f) that drives the rotation can have fluctuations, in 
which case the rotation angle varies across the sky [3H6]. 
Furthermore, if <p is some massless scalar, not necessarily 
quintessence, its background value does not necessarily 
evolve, and the uniform component of the rotation an- 
gle may vanish. The only way to look for CB in this 
scenario is through its direction dependence. Addition- 
ally, if a(n) is measured with high significance, the exact 
shape of its power spectrum provides a window into the 



detailed physics of the new cosmic scalar (j). Currently, 
the strongest limit on a direction-dependent CB angle 
comes from AGN [14] , which constrain the root- variance 
of the rotation angle to be < 3.7°. 

In previous studies [I5j [16], a formalism was devel- 
oped to search for anisotropic CB rotation with the CMB. 
The sensitivity of WMAP data to this anisotropic rota- 
tion is expected to be competitive with that from AGN 
[T5UT7] . However, the CMB also allows individual mul- 
tipoles Cf^ to be probed — the AGN data only constrain 
the variance — and is sensitive to higher L than AGN. 
The CMB also probes CB to a larger lookback time than 
AGN. 

Here we apply the formalism developed earlier to the 
WMAP 7-year data. Within experimental precision, we 
report a non-detection of a direction-dependent cosmo- 
logical birefringence. We obtain an upper limit on all 
the rotation-angle power-spectrum multipoles C^ a up to 
L = 512. This result implies a 68% confidence-level up- 
per limit on the quadrupole of a scale-invariant power 
spectrum of \/ 'C^ / / (47r) < I , 1 . As a check, we also find 
a constraint on the uniform rotation that agrees with the 
results of Ref. [13] . 

The rest of this paper is organized as follows. In ^H[ we 
review the physical mechanism for CB. In §III[ we revisit 
the full-sky formalism to search for direction-dependent 
rotation, and discuss its implementation. WMAP data 
selection, our simulations, and the tests of the analysis 
method are described in j flV| Results are reported in 
and we conclude in §VI| Appendix [A] contains a detailed 



Here, the power spectrum is defined in the usual way, Cf^ = 

X] a LM&LM /(^L + -0' wnere a spherical-harmonic decomposi- 

M 

tion of the rotation field provides the rotation-angle multipoles, 
«LM = fYl M (n)a(n)dn. 



2 



explanation of the procedure we used to obtain an up- 
per limit of the root-mean-squared rotation angle from 
the measurement of the TE correlation in the data; Ap- 
pendix [B] contains a discussion of the geometrical proper- 
ties of the rotation-angle estimator; Appendix [C] details 
the calculation of the L-dependence of the fractional cor- 
rection for a scale-invariant power spectrum recovered 
from cut-sky maps; and Appendix [D] displays the analy- 
sis masks we used in this work. 



II. PHYSICAL ORIGIN OF COSMOLOGICAL 
BIREFRINGENCE 

Theories with a weakly broken global U(l) symme- 
try provide a natural mechanism for producing a shallow 
potential for the pseudo-Nambu-Goldstone-boson field (j). 
From a cosmological perspective, a PNGB field with this 
property is a natural candidate for quintessence, since it 
can drive epochs of accelerated expansion [3]. In addi- 
tion, many other extensions of the Standard Model of 
particle physics and ACDM cosmology abound in scalar 
fields descending from theories with shift symmetry and 
resembling the PNGB. Such fields generically couple to 
photons through the Chern-Simons term F^ v F^ v of elec- 
tromagnetism, while the underlying shift symmetry sup- 
presses all other leading-order couplings to Standard- 
Model particles [3]. This way, the existence of a new 
degree of freedom <\> could evade detection in colliders 
and other lab experiments, but could still be manifest in 
cosmology through CB. We now review in more detail 
the physical mechanism that gives rise to CB. 

The Chern- Simons-modified electromagnetic La- 
grangian reads 

2 = - l - F ^F,„ - JLjF^F^, (1) 

where F^ v is the electromagnetic field-strength tensor, 
F^ v is its dual, f3 is a coupling constant, and the mass 
M is a vacuum expectation value of the spontaneously 
broken symmetry. The dispersion relation following from 
this modified Lagrangian has different solutions for the 
left- and right-handed photon polarizations, the net effect 
being the rotation of the linearly polarized electromag- 
netic wave that propagates through the vacuum with the 
evolving field <j). The direction of polarization is rotated 
by an amount 



that depends on the total change A<p along the photon's 
path. Since M can be arbitrarly large, perhaps on the 
order of the Planck mass, the accumulated change in (j) 
must also be large in order for this angle to be measur- 
able. This motivates the use of cosmological probes in 
search for CB. There are models in which A(f> is uniform 
across the sky (giving rise only to a uniform rotation an- 
gle), as well as models in which it has anisotropics [UJB]. 



In this work, we do not focus on any particular physi- 
cal model for CB, but rather derive a model-independent 
constraint on the rotation-angle power spectrum C£ a . 



III. FULL-SKY FORMALISM AND ITS 
IMPLEMENTATION 

In this Section we review the full-sky-estimator formal- 
ism of Ref. [16] for measuring direction-dependent CB. In 
order to apply this formalism to the cut sky (after mask- 
ing out the Galaxy), all measured power spectra need to 
be corrected by a factor of ~ l// s ky 2 Unless otherwise 
noted, / s ky is calculated as the fraction of the pixels that 
the mask admits; we include this factor when appropriate 
in the following derivation. We rewrite all the relevant 
formulas in a position-space form which is numerically 
efficient for the analysis of data. 

In the presence of birefringence, the polarization field 
acquires a phase factor, 

p(n) = [Q + iU](n)=p(n)e 2ia ^\ (3) 

where Q and U are the Stokes parameters for linear po- 
larization, and tilde denotes the polarization in the ab- 
sence of birefringence, which we refer to as the "primary 
polarization" . To obtain an estimate of the phase factor 
e 2ia{n) f rom polarization field in Eq. (|3|, we require a 
tracer of pip). The primary polarization is generated by 
Thomson scattering of the local temperature quadrupole, 
so the observed temperature field T(n) may be used for 
this purpose. Due to projection effects, the local temper- 
ature quadrupoles at last scattering appear on the sky as 
a curvature of the temperature field. The estimator for 
the rotation angle then involves projecting the tempera- 
ture field into a map as a spin-2 quantity (which evalu- 
ates the curvature), and looking for correlation with the 
polarization field which varies as a function of the posi- 
tion on the sky. We review a rigorous derivation of the 
estimator in the following subsections. 



A. Rotation-induced B modes 

On the full sky, the polarization field can be decom- 
posed in terms of spin-2 spherical harmonics 2^ m (n) as 

p(n) = -J2 ( E im + iB lm ) 2 Y lm (n), (4) 

Im 



2 When multipole coefficients are calculated from a map where a 
fraction 1 — f s ^ y of the pixels is masked (i.e. signal set to zero), 
the usual full-sky expression for their variance (i.e. the power 
spectrum; see Ref. [20] or Eq. ( |13| where f s ^ y = 1) is underes- 
timated by a factor of f s ky, because the variance corresponding 
to the masked pixels is effectively zero. 



3 



where E and B modes represent polarization patterns of 
opposite parity [20], [21] . 

The primary E-m.ode polarization signal Ei m (sourced 
by the dominant scalar perturbations) is detected with 
high significance in WMAP-7 data [13], although primary 
B modes (sourced by the subdominant tensor perturba- 
tions) have only been constrained with upper limits. For 
this reason, most of the constraining power for CB in 
WMAP comes from the search for a CB-induced rota- 
tion of the primary E mode into an observed B mode. 
The induced B mode is given as [I5j [16] 



Im 



dS[p(a) e 2i «( a ) 2 F4(S) 



- p(n)*e 



-2ia(n) 



■2Y; m (n)]. (5) 



This B mode is correlated with the primary E mode 
(from which it originated), and through it also with 
the temperature anisotropics. The presence of rotation 
therefore gives rise to anomalous EB and TB correla- 
tions, and both these power spectra can be used to search 
for CB. It is, however, worth keeping in mind that indi- 
vidual multipoles of the £?-mode polarization signal are 
still noise dominated, whereas the temperature is mea- 
sured at S/N > 1, for a large number of multipoles, in 
every frequency band of the WMAP 7-year data. There- 
fore, at WMAP noise levels the temperature field makes 
a better tracer of the primary £?-mode than the observed 
E mode itself. For this reason, on most angular scales, 
the search for a TB correlation, which we focus on in this 
work, provides the best constraint on CB [16 . Assuming 
the primary polarization is a pure E mode at the sur- 
face of last scatter, the CB-induced TB correlation reads 

nans] 



(Bi m T*,, 



I'm' 



dnC'Z E 



-sin(2a)[ 2 Y Fm , 2 Y 1 * 



2 Y r 



- 2 y; 



-cos(2a)[ 2 Y Fm , 2 Yf m - _ 2 Y Fm ,_ 2 Yf m ] 



, (6) 



where we suppress the n dependence for clarity. The 
two contributions to the correlation, sin and cos, have 
opposite parities, where only terms that satisfy Z + Z' + 
L =even, and Z + Z' + L =odd, respectively, contribute to 
the sum. Power spectrum Cf E is the correlation between 
the temperature and the primary E mode, which may be 
calculated with CAMB [22]. 

So far, we have not assumed anything about the mag- 
nitude of the rotation per pixel in CMB maps. Ob- 
servations of quasars suggest an upper bound on the 
root-mean-squared (RMS) of a(n) of just a few degrees 
[14], while the measurement of the TE correlation from 
WMAP-7 data implies a somewhat weaker constraint: 
(a(n) 2 ) < 9.5° (see Appendix [a1 for details). Moti- 
vated by these results, in the rest of this paper we adopt 



a small-rotation- angle limit. The numerical results we 
present in §V A| do not depend on the validity of this 
assumption, but their interpretation as an upper limit 
of the rot at ion- angle autocorrelation Cf^ does; this sub- 
tlety is discussed in more detail in j fTV] j fV] and Appendix 

El 

In the limit of small rotation angle, only the sine term 
contributes to the observed TB which then reads 



(BlmTi*, ri 



dnC^ E a(n) 



x [ 2 ^ m ,(n) 2 r; m (n) + _ 2 y^(n)_ 2 y4(n)], (7) 

A TB correlation generated by weak gravitational lensing 
of the CMB is of opposite parity, with Z+Z'+L =odd, and 
does not represent a source of bias for measuring a small- 
rotation signal. In addition, the effect of lensing is not 
internally observable at WMAP noise levels, even with 
an optimal estimator [23]. We therefore do not consider 
lensing further. 



B. Minimum- variance quadratic estimator: ollm 

From Eq. Q, it is evident that scale-dependent bire- 
fringence induces correlations between temperature and 
polarization modes at different wavenumbers Z, V; i.e., it 
produces a statistically anisotropic imprint on the covari- 
ance matrix of the observed CMB. Each IV pair measured 
in the maps may therefore be used as an estimate of the 
rot at ion- angle multipole ollm, provided that it satisfies 
the usual triangle inequalities, L < l+V and L > |Z— Z'|, as 
well as the parity condition Z+Z / +L=even. The prescrip- 
tion for combining all IV estimates in order to produce 
a minimum- variance quadratic estimator is explained in 
detail in Ref. [16 . Here, we only present the final expres- 
sions for the TB estimator, 



« LM = n l J dnY LM (n) 



Im 



I'm' 

-f (complex conjugate) 



(8) 



where Nl is an L-dependent normalization and the 
barred quantities represent inverse-variance filtered mul- 
tipoles. For full-sky coverage and homogenous noise in 
pixel space, the expressions for these quantities read 



R - B i 

-Dim = 



C BB> 



Tl'm' = 



T,f 



(9) 



The unbarred B and T are the observed temperature 
and polarization multipoles corrected for the combined 
instrumental beam and pixelization transfer function W\ , 
and the TT and BB power spectra are analytic estimates 
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of the total signal + noise power spectrum in a given 
frequency band, 



cV 



riTT, noise /tt/2 



C, 



BB 



c, 



BB 



a 



BB, noise /jjt2 



7W- 



(10) 



In the idealized case of full-sky coverage and homoge- 
neous instrumental noise, the estimator normalization 
Nl is calculable analytically and is equal to the inverse 
of the estimator variance, 



where 



E 



(2l + l)(2P + l) (CP) 2 i )2 



4-7T 



qBBqTT ' 



(11) 



C. Power-spectrum estimator: 

Once the rot at ion- angle multipoles are measured, their 
autocorrelation can be estimated as 



1 



/ sky (2L + l) 



M 



(13) 



This represents a sum over the (TBTB) trispectrum, 
where "T" and "B" denote the temperature and B-mode 
multipole moments. This estimator for Cf^ is non-zero 
even in the absence of CB-induced rotation, due to the 
presence of Wick contractions from the p rimary CMB 
and the instrumental noise (discussed in £111 B). They 
produce the "noise bias" C^' n ° lse 
tracted from the measurement of Cf^ 



and must be sub- 
in order to recover 
an estimate of the CB-induced signal 



l L V 
-2 2 



1 L V 

2 0-2 



(12) 



The objects in parentheses are Wigner-3j symbols. 

In the non-idealized case of real data, the simple 
inverse-variance filters (IVFs) presented above are sub- 
optimal, in the sense that the associated estimator vari- 
ance is not truly minimized. To obtain a true minimum- 
variance estimate, computationally more involved filters 
are required [23]. In practice, however, we find that the 
full-sky expressions for the estimator in Eq. ([8| provide a 
very good approximation to its behavior on the cut sky. 
Namely, the analytic expression for its variance, given 
by Eq. (11), is consistent with the full variance recovered 



from a suite o f Monte Carlo simulations (described in de- 
tail in <:IVB) when the simple IVFs of Eq. (10) are used 



in the presence of sky cuts; the appropriate correction for 
the fraction of the sky admitted by the analysis masks 
must be included in this case. This result motivates us 
to continue using the simple IVFs and the corresponding 
analytic expressions for the estimator normalization. 

The insensitivity to the presence of the galaxy masks 
that we observe here can be interpreted as a consequence 
of the following properties. First, the estimator of Eq. ^ 
is a product of inverse- variance filtered T and B maps in 
real space, which are local functions of the data. The 
IVFs are local in pixel space (they resemble Gaussians 
with a width of a few arcmins, corresponding to the 
resolution in a given frequency band), and so the mask 
boundaries remain localized after filtering. Additionally, 
the estimator is an even function of the temperature map 
(see Eq. Q — it contains a second derivative of the tem- 
perature field performed by 2^'m')> ano ^ so ^ * s re ^ a_ 
tively insensitive to the discontinuities introduced by the 
analysis mask. These properties put the rotation esti- 
mator Slm in sharp contrast with the estimators for the 
gravitational-lensing potential, where the dependence on 
the gradient of the temperature field renders the lensing 
reconstruction very sensitive to sky cuts [24] . 



aa, noise 



(14) 



For Gaussian CMB+noise fluctuations, the noise bias can 
be identified with the three disconnected Wick contrac- 
tions of the trispectrum which Cf^ probes: 



(a) TBTB 

|-FTI 

(b) TBTB 

i-Fh 

(c) TBTB. 



(15) 



where we neglect contraction (a), which only couples to 
the L = mode of Slm, and also contraction (c), as 
it is negligible 3 . In the absence of statistical anisotropy 
(i.e. for full-sky coverage and homogeneous instrumen- 
tal noise), the contraction (b) between two real fields 
with multipoles Zra, I'm' carries a set of delta functions 
Sii'timm' ? an d the realization-dependent noise bias may be 
written explicitly in terms of the observed power spectra. 
If Cf^ is evaluated by cross-correlating the /i, /2, ^3, and 
/4 frequency-band maps, the analytic expression for this 
"isotropic bias" follows from Eq. ([8]), 



/l/2 ry/3 /4 



noise, iso 



( a LM a LM )Gauss,iso 



(2Z + 1)(2Z' 

w 



47T 



(16) 



Ctr 



TTJ 1 f 3 maps qB B , f 2 / 4 maps 



(Pi 



TT, ft h qB B,f 2 f 2 C TT, f 3 / 3 qB B,Uf 4 y 



where the power spectra in the denominator of Eq. (16) 
are the simple analytic IVFs. The £7j T '/i/3maps 

^BS,/ 2 /4maps 



and 

are measured by cross-correlating data 
maps in the frequency bands j\ and fy, or f 2 and 



3 In our simulations, we verify that this term has indeed a negli- 
gable numerical contribution. 
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respectively, and corrected by the factor of l//gk y and 
-'-//sky' corresponding to the temperature and polariza- 
tion analysis mask, respectively. Most of the power in 
temperature comes from CMB fluctuations, and the B- 
mode power is mostly noise if ji — and negligible 
otherwise. Therefore, since the instrumental noise is in- 
dependent for different frequency bands, the largest con- 
tribution to the noise bias can be eliminated by cross- 
correlating estimates of Slm obtained from two different 
bands. 

In reality, we work with a masked sky which has been 
observed with inhomogeneous noise levels, and Eq. (16) 



does not provide a perfect description of the noise bias, 
although it is an excellent first approximation. This leads 
us to adopt a two-stage debiasing procedure in which we 
subtract both the isotropic bias of Eq. (16), and an addi- 
itonal Monte- Carlo-based correction, in order to correct 
for the effects of sky cuts and inhomogeneity of the in- 
strumental noise. The total noise bias C" a ' noise is the 
sum of the two contributions, 



C 



aa, noise 



c 



oca, noise, iso 



c 



aa, noise, MC 



(17) 



We estimate 



noise, MC 



from a set of WMAP realiza- 
tions generated with no birefringence signal (described 
in ^IVB), analyzed in the same way as the data it- 



self. For each realization, we calculate the appropriate 



C 



aa, noise, iso 



, • r^aa, noise, MC 

tarn C L as 



and averaging over many realizations ob- 



C 



aa, noise, MC 



/ fiaa s~v 



aa, noise, iso. 



(18) 



This two-stage procedure reduces the sensitivity of our 
estimator to uncertainties in the CMB and instrumental- 
noise model, as compared to the case where the entire 
bias is recovered from Monte Carlo analysis. With the 
two-stage procedure, the largest (isotropic) contribution 
to the bias is evaluated directly from the power spectra of 
the observed maps, and is specific to the CMB realization 
at hand; subtracting it from the bispectrum naturally 
takes care of any noise (bias) contribution that might 
arise from the uncertainty in the background cosmology 
or in the noise description used to generate Monte Carlo 
simulations. 

As we show in we find consistency of results ob- 
tained with either (i) the calculation of the trispectrum 
as a four-point autocorrelation of the maps in the same 
band, (ii) the calculation of the trispectrum from cross- 
band correlations, which have an almost negligible noise 
bias. 



IV. DATA AND SIMULATIONS 
A. Data 

Our main results are based on the full-resolution (cor- 
responding to HEALPix resolution of N^ e = 512) co- 



added seven-year sky maps that contain foreground- 
reduced measurements of the Stokes /, Q, and U param- 
eters in three frequency bands: Q (41 GHz), V (61 GHz), 
and W (94 GHz), available at the LAMBDA website [25] . 
A summary of the instrumental parameters most rele- 
vant to this analysis is provided in Table [I] We apply 



Band 


FWHM 


At [/iKarcmin] 


Ap [/iKarcmin] 


Q (41 GHz) 


34' 


316 


544 


V (61 GHz) 


24' 


387 


589 


W (94 GHz) 


22' 


467 


693 



Table I: Relevant instrumental parameters: beam full- width 
half maximum (FWHM) and approximate map noise for tem- 
perature and polarization for the three frequency bands we 
use in the analysis 25 . 



the seven-year temperature KQS5y7 mask with 78.27% 
of the sky admitted, and a polarization P06 mask with 
73.28% of the sky admitted. These masks are constructed 
to remove diffuse emission based on the data in K and Q 
bands, and on a model of thermal dust emission. Point 
sources are masked based on a combination of external 
catalog data and WMAP-detected sources. (For more 
information about the exclusion masks, see Appendix [D] 
and Ref. [27].) 



B. Simulations 

We produce a suite of simulated WMAP observations, 
both to test the normalization of our o^m estimates as 
well as to estimate their variance for the subtraction of 
C 



aa, noise 



in Eq. (14). We produce simple simulations of 



the WMAP data with the following procedure: 

1. Generate CMB-sky temperature and polarization 
realizations for the best-fit "LCDM+SZ+ALL" 
WMAP-7 cosmology of Ref. [25]. 

2. Convolve the simulated CMB skies with a symmet- 
ric experimental beam. For the WMAP band maps 
we use an effective beam calculated as the average 
beam transfer function for all differencing assem- 
blies at the given frequency. 

3. Add simulated noise realizations based on the pub- 
lished /, Q, U covariance matrices within each pixel. 
We do not make any attempt to generate noise with 
pixel-to-pixel noise correlations, although we do ex- 
clude multipoles with I < 100 from our analysis, 
as this is where most of this correlated noise re- 



sides. In jVA we demonstrate the consistency of 



C^ a estimates constructed from auto- and cross- 
correlations of maps with independent noise real- 
izations, and so are justified in neglecting correlated 
noise in our analysis. 
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We do not include Galactic foreground residuals or unre- 
solved point sources in our simulations, but we address 



their possible impact on our results in jVC 



C. Test runs 

In order to demonstrate the recovery of the CB signal 
using the minimum- variance estimator formalism and the 
de-biasing method discussed in the previous sections, we 
generate a suite of simulations that include a CB signal, 
i.e. where the polarization maps are rotated by realiza- 
tions of a scale-invariant power spectrum of a, 



QOLOL 



AC 



aa, fiducial 



A 



131deg 2 
L(L+1)' 



(19) 



where we choose the amplitude of this fiducial model so 
that it gives a S/N ratio of order 1 at low L for WMAP V 
band, and an RMS rot at ion- angle on the sky of 10°, satis- 
fying the small- angle approximation 4 . We apply analysis 
masks to each simulated map, and then analyze the map 
cross-correlations (as described in previous sections), re- 
covering Slm multipoles; we then compute the power 
spectrum using Eq. 13 Due to the interaction of the 



power distribution at different scales in the map with 
the geometry of the analysis mask, the / s k y factor is in 
principle a function of the multipole moment L, which 
typically starts smaller than the average 5 value at low 
L's, and converges towards the average value at high L's. 
Since most of the signal for this model (which we come 
back to in ^VB) comes from low L's, we evaluate the ex- 



act L dependence, and substitute the f s ^ y (L) function in 
Eq. 13 (for more details on / s k y (L), see Appendix |C|) . In 
Fig /|Tf show the results of this test, comparing the input 
Qaa power spectrum to the mean of the reconstructed 
power from a large number of simulations, and demon- 
strate a successful recovery of the signal. In the following 
Section, we apply the same signal-reconstruction method 
to WMAP-7 data. 

To conclude this section, we note one subtlety nec- 
essary for the correct interpretation of the results of our 
analysis. The expression for the estimator of Eq. ^ only 
recovers the rotation-angle multipole in the small-angle 
regime. In the general case of arbitrarily large rotation, 
Eq. ([8| provides an exact estimate of the multipoles of an- 
other "observable" quantity: |sin[2a(n)]. Strictly speak- 
ing, our de-biasing procedure also relies on the small- 

i j • s^aa, noise, MC • i i , i 

angle approximation, since C L is calculated 



4 It is important that the model satisfies the small-angle approxi- 
mation, as our calculation of the bias from Monte Carlo analysis 
is based on the null-assumption. In regime where this approxi- 
mation is not satisfied, higher-order corrections will be necessary 
to recover the rotation-anlge power spectrum from the measured 
< TBTB > trispectrum. 

5 The "average" here is the usual / s k y fraction associated with a 
mask, equal to the fraction of the pixels that the mask admits. 



from a suite of null simulations. It is therefore necessary 
to inquire which regime corresponds to a particular model 
of rotation before interpreting our results as a constraint 
on such a model. However, the fiducial model we use as 



an example here (and which we come back to in £VB) 
satisfies this assumption (producing an RMS rotation of 
~ 10°). In this particual case, the difference between the 
two power spectra of a and of \ sin [2a] is mainly con- 
tained in the 15% difference in their amplitudes. It thus 
is possible to recover the rot at ion- angle power spectrum 
by simply rescaling the measured power — the fact we use 



in i V B to constrain this model from WMAP data. 



0.0201 



0.015 



— - input a 

— input lsin(2a) 

— measured a 



, 0.010 



0.005 



0.000 




Figure 1: The recovery of the CB signal with the analysis 
method presented in §III| is demonstrated using a suite of sim- 
ulations that include realization of a rotated sky. Blue dashed 
line is the input signal power-spectrum of C£ a , red thick line 
represents the power spectrum of |sin[2a(n)], and the thin 
black line is the mean recovered power from the simulations; 
the gray region is a lcr-confidence-level interval calculated 
from the same suite of simulations. 



V. RESULTS 

A. Model-independent constraints 

Before continuing, let us first clarify our notation. The 
rotation-angle power spectra are marked with four fre- 
quency bands as [/1/2] L/3./4] • This means that the two 
estimates of ollm needed to evaluate the power spectrum 
are obtained by cross-correlating band f\ with / 2 , and 
fs with /4, respectively. Here, the temperature multi- 
poles are measured from fi and ^3, and the B modes 
are obtained from the maps in fi and bands. We 
measure five different band cross-correlations: [VV] [VV] , 
[QV][QV], [QQ][W], [WV][WV], and [WW][VV], but 
since the results for all of them are qualitatively the same, 
here we only show plots for a characteristic subset. 
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Figs.[2j[3j and[4]show the measurement of the rotation- 
angle autocorrelation, before and after debiasing, and dif- 



ferent components of the noise bias described in £ 111 C 
The blue and gray areas in the middle panels represent 
la and 3a confidence-level intervals, respectively, derived 
from the null-hypothesis (no rotation) Monte Carlo anal- 
ysis described in §IVB| We see no significant deviations 
from zero in any of the five band cross-correlations — our 
results are consistent with ollm = to within 3<j, at each 
multipole in the range from L = 0toL = 512. We bin 
the power and list the measurements for all multipoles 
in Table [TTJ As an additional consistency check, the up- 
per limit we obtain on the uniform rotation angle, given 
as a = aoo/V&r, is in good agreement with previous 



WMAP results [13]; see Table III 



As we pointed out in §IVC[ in the general case of ar- 
bitrarily large rotation, our method provides an exact 
estimate of the autocorrelation of the following quantity: 
|sin[2a(n)], rather than the rotation angle itself; when 
the small-angle approximation is satisfied, this quantity 
and its power spectrum assymptote to a(n) and re- 
spectively. In order to establish the regime corresponding 
to a particular model, we note that the RMS fluctuation 
typical of realization of a power spectrum C^ a is given 

by 



<a(n) a ) 1/a 




(20) 



In the event of a breakdown in the small-angle approx- 
imation, the values in Table [TT] should be interpreted as 
constraints on the autocorrelation of |sin[2a(n)], rather 



than a itself. Evaluating Eq. ( 20 ) for the uncertainty lev- 
els quoted in Table [TT] would erroneously lead to a conclu- 
sion that a large RMS rotation is allowed by the WMAP 
data; we show that the upper limit on the RMS rotation 
is roughly 9.5° (see Appendix [A]), and we again note that 
previous studies of quasar data imply an even stronger 
constraint of ~ 4° [13] . 



B. Constraints on a scale- invariant power spectrum 



The null result shown in £ V A can be translated into an 



upper limit on the amplitude of any model of rotation. As 
a generic example, we focus on a scale-invariant rotation- 
angle power spectrum of Eq. ( 19 ). The best-fit amplitude 



is evaluated from all multipoles in the range < L < 512, 
using a minimum-variance estimator [26] 



/^aa, fiducial Piaa 

® A 



L 



var(C£ a ) 



(21) 



where 



O A 



(C 



tt a, fiducial ^2 N 



-1/2 



var(C£ a ) 



(22) 



L bin 


[W][W] 


[QV][QV] 


[QQ][W] 


[WV][WV] 


[WW][W] 


26 


2.65±1.87 


1.61±2.44 


1.05±1.62 


0.72±2.03 


-0.43=bl.34 


77 


1.86±2.58 


0.70zb2.84 


1.57±2.36 


1.03±2.70 


0.17±2.04 


128 


1.07±1.33 


1.00±1.36 


0.27zbl.l7 


3.04±1.35 


0.96±1.02 


179 


1.40±1.49 


-l.29zbl.65 


-0.31zbl.15 


-0.40dbl.48 


0.66±1.13 


230 


-l.90zbl.76 


-4.47zbl.96 


1.87±1.36 


-3.36±1.97 


-0.69±1.33 


282 


4.31±2.23 


3.17zb2.42 


2.04±2.21 


2.14±2.42 


-0.20zbl.90 


333 


1.98±2.39 


-0.25zb2.60 


4.59=bl.80 


2.62±2.45 


-l.llzbl.96 


384 


0.81±1.78 


-1.71±1.93 


1.97zbl.51 


1.22zbl.71 


1.93zbl.52 


435 


-0.40±1.64 


-0.19=bl.80 


1.53zbl.26 


-1.03±1.74 


-l.65zbl.30 


486 


3.22±1.75 


0.78±1.93 


1.02±1.39 


2.69=bl.84 


-0.28±1.27 



Table II: Results for the measurement of Cg" [degrees 2 ] are 
listed, as recovered from five different band cross-correlations. 
The la confidence intervals are calculated with a suite of 
Gaussian sky simulations, described in §IVB[ The results are 
binned, with the central L value of each bin listed in the table; 
the width of each bin is AL ~ 51. 



[/1/2] 


a±la [°] 


[VV] 


-0.9 zb 2.3 


[QV] 


-0.5 zb 2.4 


[QQ] 


0.9 zb 2.8 


[WV] 


-2.2 zb 2.4 


[WW] 


-1.8 zb 2.7 



Table III: Uniform-rotation angle a with a la confidence in- 
terval, from five cross-band correlations of WMAP-7; the cor- 
rection factor of l//sky is applied to each measurement here. 
The uncertainties are consistent with the zbl.4° uncertainty 
on the uniform-rotation angle reported by the WMAP team 
[T3] for a joint analysis of the Q, V and W-band data, after 
accounting for the fact that we analyze the bands individually 
(resulting in slightly larger error bars). 



is the analytic expression for the variance of A, and 
var((7£ a ) is the variance of the null- hypothesis rotation- 
angle power spectrum, estimated from a suite of null- 
hypothesis simulations. We note that the measured 



(jaxx h ave been corrected by f s ^ y (L) (see <TV C and Ap- 
pendix [Cj the correction is calculated specifically for this 
model) only in this Subsection; for the presentation of 
the model-independent results, we use the average value, 
/sky ~ 0.68. Most of the constraint here comes from low 
L's; 50% of the sum in Eq. (22) comes from L = 1, and 
90% from L < 10. 

Even though the analytic expression above provides a 
good estimate of the statistical variance, because the con- 
straint comes primarily from low-L modes the probability 
distribution of A is significantly no n- Gaussian. To cap- 
ture this non-Gaussianity in our analysis, we again gen- 
erate a suite of null-hypothesis Monte Carlo simulations 
and recover the 68% and 99% confidence-level intervals 
from these simulations. The corresponding probability 
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Figure 2: Top panel: Measurement of the rotation- angle 
power spectrum from V band, shown before debiasing, along 
with the components of the noise bias: Monte-Carlo mea- 
surement of the null- hypothesis mean (C^ a ) (solid green), 
isotropic noise bias (blue dashed) , and the mean isotropic bias 
(magenta dashed). Middle panel: The same power spectrum 
after debiasing, with la and 3a confidence interval. Bottom 
panel: binned version of the middle-panel power spectrum. 
The results are consistent with zero within 3a. 

distributions for A are shown in Figure |ij 

The best-fit values for the quadrupole amplitude Cg 01 
and associated confidence intervals are listed in Table 
|IV| consistency with zero is apparent within 3a for all 
band-cross correlations we analyzed. The tightest con- 
straint on the quadrupole amplitude of a scale- invariant 
rot at ion- angle power spectrum comes from [WW] [VV] ; it 
is VQTA 4 ^) < 1° with 68% confidence 6 . 



Figure 3: Same as Figure^ for [f 1 f 2 ] [/ 3 / 4 ] = [WW] [VV]. 



[/1M/3/4] 


CF ± la(±3a) [deg 2 ] 


[VV][VV] 


/, + 15.8/+79.0x 
11 -^-16.9 1-27. 77 


[QV][QV] 


9Q £+18.8/+70.3x 
zy,u -18 .3 V-33 A) 


[QQ][vv] 


1Q o + 14.3/+51.6\ 
J-y-8-13.9U46.eJ 


[WV][WV] 


o + 15.9/+79.0x 
lu -°-16.9 \-27.7) 


[WW][VV] 


n +14.0/+43.3\ 
°- u -13. 91-42.97 



Table IV: Measurement of the quadrupole amplitude of a 
scale- invariant rotation-angle power spectrum for different 
cross-band correlations, with 68% and 99% confidence-level 
intervals, recovered from a suite of null-hypothesis simula- 
tions. Consistency with zero within 3a is apparent for all 
band cross-correlations, and the tightest constraint comes 
from [WW][VV]. 



Note that the conversion between the amplitude A and the 



quadrupole is C2 = A x 131deg 2 /6. 
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[W V][W V] 




[W V][W V] 




Figure 4: Same as Figure^ for [/i/ 2 ][/ 3 /4]=[WV][WV]. 



C. Potential systematics 



In addition to the statistical error reported here, there 
is also a systematic error for the measurement of the uni- 
form rotation angle, owing to uncertainty in the detec- 
tor polarization angles [13]. This systematic uncertainty 
should only apply to the monopole of a. The direction- 
dependent part is only sensitive to the extent that it 
affects the statistical noise bias, and this is mitigated 
by our data-dependent debiasing procedure. There are, 
however, other sources of systematic error that can po- 
tentially bias our estimates and add uncertainties to the 
rotation- angle measurements. In this Section we investi- 
gate the impact of asymmetry of the instrumental beams, 
unresolved polarized point sources, and foreground resid- 
uals from unremoved/unmasked Galactic emission. 



1. Beam Asymmetries 

The fast spin and precession rates of the WMAP scan 
strategy, coupled with the yearly motion of the satellite 
around the ecliptic, enforce that any bias to S^m origi- 
nating from scan-strategy related systematics, like beam 
asymmetry, must be confined to M = modes in eclip- 
tic coordinates [29]. Furthermore, the smoothness of the 
scan strategy on large scales (dictated by the 85-degree 
opening angle of the detectors, and the large 22.5-degree 
amplitude of the hourly satellite precession), ensure that 
any such bias falls off quickly as a function of L. The esti- 
mates of C2, as we have discussed in the previous section, 
are most sensitive to the low-L modes of C£ a , so to test 
for the presence of beam-asymmetry contamination, we 
rotate our coordinate system to ecliptic coordinates, and 
re-derive a constraint on C2 from L < 10, M = modes. 
We see no departure from the null hypothesis in this case 
where it should be maximal, and so conclude that beam 
asymmetries are not a significant source of bias for our 
measurements. 



2. Unresolved point sources 

To test the impact of unresolved point sources on our 
results, we repeat the analysis after unmasking the por- 
tion of the maps which are associated with detected point 
sources (but not Galactic contamination). Compared to 
our fiducial analysis, the measurement points for Cf^ 
shift by < la, where a represents the statistical error 
from our foreground-free Monte Carlo analysis; see Fig- 
ure [6] This shift provides a conservative upper limit 
on the systematic uncertainty that point-source residuals 
can produce, assuming that the bright detected popula- 
tion has comparable polarization properties to those of 
the fainter sources. The unresolved point-source power 
at WMAP frequencies is dominated by unclustered radio 
sources, with fluxes close to the detection threshold, and 
so this is a reasonable assumption. We note that there 
is no overall bias, as the direction of scatter does not ap- 
pear to be correlated for different multipoles. Of course, 
the contribution of radio point sources to the map is a 
steep function of the flux cut, and by unmasking all de- 
tected point sources our estimate of potential bias and 
uncertainty is overly conservative. Given a model for 
radio-source number counts dN/dS, we can scale these 
results to the levels of contamination expected at the 
actual WMAP source detection threshold of (conserva- 
tively) ~ 1 Jy. Any bias AC£ a (which we do not see 
evidence for, even in the unmasked map) will scale with 
the point-source trispectrum as 



AC aa 



pS cu t 

Js=o 



(23) 
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WWW: Histogram of best-fit amplitude? 




WVWV: Histogram of best-fit amplitudes 




VWV: Histogram of best— fit amplitudes 




Figure 5: Probability distributions of the best-fit amplitude A of the scale- invariant rotation power C£ a recovered from a 
suite of null- hypothesis realizations are shown for some of the cross-band correlations. The gray-shaded areas denote a 68% 
confidence interval around a mean value; the red vertical line represents the measurement of the best-fit A for a given band-cross 
correlation. We find consistency with zero within 3a for all five measurements. 



while the uncertainty on Cf^ will scale with the point- 
source power as 

'(<*■)«* < 24 > 

Evaluating these terms for the dN / dS model of Ref. [28] , 
we find that A and a should be suppressed by factors 
of 0.005 and 0.06 respectively when moving from a flux 
cut of 10 Jy (no masking) to 1 Jy. We find even smaller 
(though comparable in magnitude) results using the sim- 
pler dN/dS oc 5'- 215 model of Ref. [30,. This implies 
that any bias from unresolved sources should be com- 
pletely negligible, and any increase in uncertainty due 
to their contribution to the observed power should be 
< 0.1a, where a represents the statistical error from 
our point-source-free Monte Carlo analysis. In conclu- 
sion, we expect the unresolved point sources to produce 
a negligible systematic uncertainty in the measurement 
of Cf\ 



3. Foreground residuals 

An additional conceivable source of systematic uncer- 
tainty might result from Galactic foregrounds. To ex- 
plore the extent to which such uncertainty might affect 
our results, we perform two tests. In the first, we repeat 
our analysis on non-foreground-reduced maps, to test 
the effect of the presence of unsubtracted foregrounds. 
In the second, we repeat the foreground-reduced analy- 
sis using a mask which excludes a larger fraction of the 
low-Galactic-latitude sky. We construct this conservative 
mask by combining the fiducial KQ85y7 analysis mask 
with the extended mask of Ref. [25 , and additionally 
masking out pixels with Galactic latitudes in the range 
of ±40°; the mask admits only about 33% of the sky, ap- 
proximately half the sky admitted in our fiducial analysis 
(where ~ 68% of the sky is admitted; see Appendix |D|) . 
The function of this test is to explore the effect of resid- 
uals left by the foreground subtraction procedure, which 
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1 \ 




. • Foreground reduced, analysis mask 

O Foreground reduced, point sources unmasked 



100 200 300 400 500 



L 

Figure 6: Measurement of from [WW][VV]. Results 

shown in black are obtained by using the analysis mask that 
covers all the point sources brighter than ~ Uy, while the re- 
sults in red (empty circles) are obtained after unmasking all 
the point sources. There is no apparent bias and the difference 
in every bin is less than the statistical uncertainty, despite the 
extreme variation in the source contamination. Scaling argu- 
ments in §VC 2| imply that the unresolved point sources have 
a negligible contribution to the estimated measurement un- 
certainty for the most constraining band-cross correlation. 



should be stronger close to the Galactic plane. The mea- 
surement of C^ a is scaled appropriately to account for 
the fractional sky coverage and the results from the two 
modified analyses are compared with the results of the 
fiducial analysis in Figure [7] In the first case, the change 
in the measurements is negligible compared to the sta- 
tistical uncertainty. In the second case, the scatter be- 
tween the two results is consistent with the difference in 
sky coverage (producing up to 40% larger scatter for the 
extended- mask data points). The measurements show 
no apparent bias in either case. These results imply that 
foregrounds and foreground residuals are not likely to 
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make a large systematic contribution to our estimated 
statistical uncertainty, at least for the case of the most 
constraining band-cross correlation [WW] [VV] . 
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Figure 7: Measurement of Cg" from [WW][W] band cross- 
correlation, with the corresponding statistical uncertainty ob- 
tained from a suite of null-hypothesis foreground-free simula- 
tions. Black filled circles are measurements obtained from 
the foreground-reduced maps after applying the fiducial anal- 
ysis mask, and they represent our fiducial results of Figure 
[3] The colored data points (and the associated error bars) 
correspond to the two test analyses: the green x's are ob- 
tained from the maps prior to foreground subtraction, but 
using the fiducial mask, while the red empty circles are mea- 
surements obtained from foreground-reduced maps after ap- 
plying an extended mask (with an additional ~ 35% of the sky 
covered around the Galactic plane) . The results from the two 
tests show no apparent bias. For the case of foreground-non- 
reduced analysis, the difference from the fiducial measure- 
ments is negligible compared to the statistical uncertainty; for 
the extended-mask case, the scatter between the two results 
is consistent with the difference in sky coverage. This implies 
that foregrounds and foreground residuals should not have a 
drastic impact on the estimated measurement uncertainty. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we implement the minimum-variance 
quadratic-estimator formalism of Ref. [16] to search 
for direction-dependent cosmological birefringence with 
WMAP 7-year data. We derive the first CMB mea- 
surement of the rotation-angle power spectrum in the 
range L = — 512, finding consistency with zero at 
each multipole, within 3<r. We estimate an upper limit 
on each power-spectrum multipole by simulating a suite 
of Gaussian sky realizations with no rotation, including 
symmetric beams and noise realizations appropriate for 
each WMAP frequency band, and also Q-U correlations 
and sky cuts. We investigate the impact of foregrounds 



and polarized diffuse point sources on the reported con- 
straints and come to the conclusion that they are not 
significant sources of systematic error for the rotation- 
angle estimates. Finally, we use the null result to get 
a 68% confidence-level upper limit of yC^7 (^0 < 1° 
on the quadrupole of a scale-independent rotation- angle 
power spectrum. Even though the CMB constraint turns 
out to be comparable to that derived from quasar mea- 
surements, the CMB analysis has a significant advantage: 
it provides a measurement of the rot at ion- angle power at 
each individual multipole L and has better sensitivity to 
models with significant power at high multipoles. 

The same formalism we use in this work can be 
applied to forecast the sensitivity of upcoming and 
future-generation CMB satellites to detecting direction- 
dependent cosmological birefringence. With 7 years 
worth of integration time with WMAP, we are able to 
constrain the uniform component of the rotation to less 
than about a degree; it will be interesting to see the re- 
sults of this analysis method applied to the upcoming 
data from Planck satellite [31], where the sensitivity to 
rotation angles on the order of a few arcminutes [16] is 
expected. 
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Appendix A: Constraints on rms rotation from 
WMAP 

If the primordial B mode is small compared to the 
primordial E mode, and the rotation field is independent 
of the CMB, the measured TE correlation reads (see also 
Figure [8} 



a 



TE 



-<cos[2a(n)])Cr, (Al) 

where the mean is taken over all realizations of the ro- 
tation field, and it does not depend on the direction n. 

In the case the probability distribution for a is a Gaus- 

1/2 

sian centered at zero and with a width (a(n) 2 ) , the 



expectation value in Eq. (Al) is simply related to the 
pixel- variance of a, 

e -2(a(n) 2 )_ 



(cos[2a(n)]) 



(A2) 



Therefore, an estimate of this expectation value and its 
uncertainty, obtained from the TE measurement as com- 
pared to the primordial power sectrum Cj E provides an 
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upper limit of the rotation- angle pixel- variance. Adopt- 
ing the expressions for a minimum-variance estimator 
and its variance (analogous to those of Eqs. (21) and 
([22])), we obtain: (cos[2a(n)]> = 0.997 ±0.050 (note that 
this constraint follows from the 21a confidence- level de- 
tection of the TE correlation reported by Ref. [13]), im- 

1 /2 

plying an upper limit on the rotation rms (a(n) 2 ) < 
9.5°. 
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Figure 8: TE correlation measured from WMAP-7 data 
(black) is compared to the primordial power spectrum, gener- 
ated using CAMB for the best-fit cosmology with no rotation 
(red, solid); both power spectra are obtained from Ref. |25| . 
The uncertainty on this measurement (gray) leaves room for 
a maximal rotation-angle rms of roughly 9.5°. 



power-spectrum kernel as well as the geometric Wigner- 
3j contributions to the CB kernels in the IV space from 
Eq. (|8|; see Figs. [9] and 10 The structure of the power- 
spectrum kernel originates from the polarization and 
temperature power spectra; the terms that correspond to 
the acoustic peaks in the primordial TT and TE power 
spectra have the largest contribution to the sum over IV. 
The geometric weight dictates the shape of the l,V,L 
triangles which are generated by CB at a scale L. The 
terms where either I or V is close in value to L have the 
largest contribution. The combination of the geometric 
weight and the power-spectra weight dictates the size of 
the noise bias at any particular L. The interplay of the 
two, for example, produces a peak at L ~ 270, apparent 
in all the plots of the noise bias presented in this work. 
The local maximum in the variance of Cf^ also appears 
at this scale. 

Appendix C: f s ^ y (L) 



In order to evaluate the exact L-dependence of / s k y 
used to reconstruct the scale-invariant rotation-angle 



power from the CMB maps (see SIVC), we generate a 
large number of a(n) realizations of the power- spectrum 
model of Eq. ( 19 ), mask the sky with the fiducial analysis 



mask, and then recover the input power spectrum in the 
usual way, i.e. take the pseudo-Cx of the masked map. 
The / S ky (L) shown in Fig. 11 is the average ratio of the 



output to the input power, as a function of the multipole 
moment L. 



Appendix D: Analysis masks 



Appendix B: Visualizing the CB kernels 

To illustrate the shape in harmonic space of the sta- 
tistical anisotropy introduced by CB, we plot here the 



Here we vizualize all the analysis masks we used in this 
paper: the fiducial temperature and polarization masks, 
as well as the two test masks used in §V C| 
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Wigner 3j kernel: Hi (2Z + 1)(2Z' +1)/(4tt) [L=200] 





Figure 9: The Wigner-3j geometric factors in the summands of Eqs. (|8| and (16) for three different values of L are shown. 



The geometric factor is non-zero only in the region of the 11' space where the triangle inequalities are satisfied. The dominant 
contribution comes from the triangles in which L ~ /, or L ~ l\ i.e. where either the temperature or the polarization mode has 
a length scale comparable to the length scale of the rotation-angle mode. 



An alytic kernel: (cfV /(C^cF) [V ba nd] 
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Figure 10: T he p ower-spectrum kernel of the summands in 
Eqs. (pi and (16) is shown. The types of triangles that con- 



tribute the most to the isotropic bias of are set by the 
geometric properties of spin-2 Wigner-3j symbols illustrated 
by the kernel shown in Figure [5] which is modulated by this 
kernel to produce summands in the expression for the bias. 



Figure 11: The L dependence of the / s k y factor used for 
the reconstruction of the scale-invariant rotation-angle power 
spectrum in £IVC and [VB The horizontal (green) line at 
/sky = 0.68 represents the fraction of pixels admitted by the 
mask. 



A. De Rosa, F. Fi nelli and F. Paci, JCAP 1202, 023 
(2012) |arXiv:1107.5548| [astro-ph.CO]]; R. Antonucci, 
|nep-ph/0110206| 
[9] A. Lue, L. M. Wang and M. Kamionkowski, Phys. 
Rev. Lett. 83 1506 (1999) [arXiv:astro-ph/98 12088 1; 
N. F. Lepora, arXiv:gr-qc/9812077| G. -C. Liu, S. Lee 
and K. -W. Ng, Phys. Rev. Lett. 97, 161303 (2006) 
|astro-ph/0606248| , 
[10] B. Nodland and J. P. Ralston, Ph ys. Rev. Lett. 78, 
3043 (1997) [arXiv:astro-ph/9704196'] ; T. J. Loredo, 
E. E. Flanagan and I. M. Wasserman, Phys. Rev. D 
56, 7507 (1997) |arXiv:astro-ph/970 6258 ; D. J. Eisen- 
stein and E. F. Bunn, Ph ys. Rev. Lett. 79, 1957 (1997) 
|arXiv:astro-ph/9704247l ; J. F. C. Wardle, R. A. Per- 
ley and M. H. Cohen, Ph ys. Rev. Lett. 79, 1801 (1997) 
|arXiv:astro-ph/9705142 ; J. P. Leahy, arXiv:astro- 
ph/9704285| S. M. Carroll and G. B. Field, Phys. 



Rev. Lett. 79, 2394 (1997) |arXiv:astro-ph /9 704263 ; 
A. Cimatti, S. di Serego Alighieri, G. B. Field and 
R. A. E. Fosbury, Astrophys. J. 422, 562 (1994); 
S. d. S. Alighieri, F. Finelli an d M. Ga laverni, Astro- 
phys. J. 715, 33 (2010) [arXiv:1003.4"823l [astro-ph.CO]]; 
P. P. Kronberg, C. C. Dyer, and H.-J. Roser, Astro- 
phys. J. 472, 115 (1996); D. Horns, L. Maccione, A. Mi- 
rizzi and M. Roncadelli, |arXiv: 1203.2184 [astro-ph.HE]; 
L. Shao and B. -Q . Ma, Phys. Rev. D 83, 127702 (2011) 
[arXiv: 1104.4438] [astro-ph.HE]]. 
[11] H. C. Chiang, P. A. R. Ade, D. Barkats, J. O. Battle, 
E. M. Bierman, J. J. Bock, C. D. Dowell and L. Duband 
et a/., Astrophys. J. 711, 1123 (2010) [arXiv: 0906.11811 
[astro-ph.CO]]. 

[12] E. Y. S. Wu et al [QUaD Collaboration], Phys. Rev. 

Lett. 102, 161302 (2009) [arXiv: 08 1 1 . 0618] [astro-ph] ] . 
[13] E. Komatsu et al. [WMAP C ollaboratio n], Astrophys. J. 

Suppl. 192, 18 (2011) ;arXiv:1001.4538| [astro-ph.CO]]. 



14 





Figure 12: Fiducial temperature-analysis mask (with ~ 78% of the sky admitted), fiducial polarization-analysis mask (~ 73%), 
extended mask (~ 33%), and the temperature-analysis mask with point sources uncovered (~ 82%). Fraction of the sky 
admitted of the combined fiducial masks for polarization and temperature is ~ 68%, which is twice the sky admitted as 
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